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Abstract 

We construct models of self-interacting scalar fields whose BPS solutions exhibit 
kink profiles which can be continuously deformed into two-kinks by varying one of the 
parameters of the self-interacting potential. The effective models are obtained from 
other models with two interacting scalar fields. The effective models are then applied 
in a brane-world scenario where we analyze the consequences of the thicker branes in 
the warped geometry and in the localization of gravity. 

1 Introduction 

Nonlinear systems, particularly those that exhibit topological solutions are very important 
for modeling many physical systems from condensed matter to high-energy physics. One 
example that has been receiving considerable attention is that of topological structures in 
multidimensional warped space-time when one considers scalar fields coupled to gravity. In 
this context, one of the key points is the possibility of localizing gravitons in thin branes 
and reproducing effectively the four- dimensional gravity [1]. In (4,1) dimensions it is shown 
that thick branes can also localize the gravity [2]. In fact, the subject of thick branes in 
the context of brane- worlds has received a considerable amount of attention [3j . Some years 
ago it was observed that some kinds of models with two interacting scalar fields in a warped 
geometry can be used to describe the splitting of thick branes due to a first-order phase 
transition [4] . A few years ago, it was shown that Bloch branes are solutions of a model with 
two interacting scalar fields which can be naturally incorporated in a supersymmetric theory 
[5] . Later, that same model was shown to exhibit branes with richer structures [6j than those 
found in the reference |5], including critical and degenerate branes. Moreover, contrary to 



* e-mail: hott@feg.unesp.br 



1 



the scenarios presented in [1] and [5] where the sphtting of the branes is controlled by a 
coupling constant presented in the interaction potential, it was shown [6] that the thickness 
of the branes is controlled by a parameter, called degeneracy parameter, which is not present 
in the Lagrangian density of the model. Instead, it is one of the constants of integration 
of the orbit differential equation which relates both fields. Recently, a model with only one 
scalar field which incorporates thicker branes was proposed by Dutra [7J. It is characterized 
by a non-polynomial interaction potential with coupling constants that control the thickness 
of the brane. This last property, besides the fact of exhibiting a non-polynomial interacting 
potential, is shared with the p-model introduced in the references [8j and [9J. In the later 
one, the changing of the parameter p in discrete jumps implies into a changing of model 
and thicker branes appears only for some values of p, whereas, in the model introduced in 
[7], the parameter that controls the deformation of the brane is a coupling constant of the 
model and as such, its variation do not modify the structure of the model. 

The purpose of this work is to construct nonlinear models, in classical field theory, with 
only one scalar field, from models with two interacting scalar fields and that exhibit solutions 
with two-kinks profiles. Two-kinks solutions yield thicker branes in brane-worlds scenarios 
or Bloch branes whose internal structures are somehow incorporated in the parameter that 
controls the thickness of the brane. The non-polynomials effective models we construct 
belong to the same class of models than that proposed by Dutra [7j, where the kink solutions 
can be continuously deformable into two-kinks by varying one of the parameters of the 
effective potential. Moreover, the same parameter that controls the thickness of the brane 
can be thought as depending on the temperature in such a way that a first-order phase 
transition characterized by the brane splitting can happen. The model with two interacting 
scalar fields we consider here is the same whose consequences in brane-worlds scenarios were 
already studied in the references [5] and [6j. In order to construct the effective models, 
we apply a general orbit equation relating both fields [10] to eliminate one of the fields in 
favor of the other. The resulting effective model is applied to the study of the localization 
of gravity by thick branes in (4,l)-dimensional warped space-time. In the next section we 
introduce the kind of models with two interacting fields we are interested in and, after that, 
we show how the orbit equation can be used to construct models with only one scalar field. 
One particular example is studied in detail. In the third section we apply one of the effective 
models in brane- world scenario and analyze the consequences of thick branes over the warp 
factor and on the localization of gravity. In the conclusions section we remark on possible 
applications of the effective models, particularly in the phenomenon of brane splitting. 

2 Models with only one scalar field constructed from 
models with two interacting scalar fields 

The models with two interacting fields we consider here are described by the Lagrangian 
density 

£ = ^d^<pd,<P + ^d^xd.X - Vi<P, x), (1) 
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with the potential written as 
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V{<P,X) = ^{W} + W^). (2) 

where and are the derivatives of some function W{(f),x), called superpotential, with 
respect to the fields (p and respectively. 

It is shown that for potentials written in terms of a superpotential, as in equation ([2]), 
the static solutions of the first-order equations 

are those that minimize the energy of the system, the BPS energy [TT], and are also solutions 
of the static equations of motions [T^. Based on the first-order differential equations one 
can realize that dx is a kind of invariant and the following, in general nonlinear equation 

furnishes a relation between the classical static solutions, called the orbit equation, which 
can be solved analytically depending on the model under consideration [TU]. For the cases 
in which the general orbit equation can be found, we can write the scalar field x terms of 
the scalar field 0, that is x = fi.'P)^ ^-^id it can be used to eliminate x terms of (j) in the 
Lagrangian density ([1]), that is 

^ = ^(l + /JW5M0-m, (5) 

where /</> = ^ = ^ = ^ with both, and W^, written in terms only of the field (j). V{(j)), 
is the potential V{(j), x) with x written in terms of (j). Due to the peculiar structure of the 
potential given in ([2]) and also to the differential equation (jll), we have 

v{4>) = \{i + fi)wl (6) 

The equation of motion for this single scalar field model is given by 

(1 + fl)d,d^<i> + UU^d,<pd'^<p + ^ = 0, (7) 

and the energy associated with the static classical solutions is expressed as 

E = I + + (8) 

where 0' stands for d(f)/dx.ln order to find the minimum energy we note that the energy can 
be rewritten as 

^ I dx [((1 + ± V2VY T 2(1 + flY/'<P'V2v] , (9) 



E = - I dx 
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and, consequently, the classical solutions with minimum energy satisfies the first-order dif- 
ferential equations 

2V 

<P' = T(t-7i)^/' = TW^, (10) 

where = dW{(j), x)/d(p with x replaced by It is easy to show that the solutions of 

the first-order differential equations are also solutions of the static equation of motion 

(1 + + uu^f = ^. (11) 



and the BPS energy is given by 

Ebps = 



\W{^, /(0))(oo) - Wi^, /(0))(-oo)|, (12) 



where the classical solutions (f){x) are to be taken at ±cxo. 

This procedure resembles the one carried out in the reference [12] to prove the equivalence 
between sine-Gordon, Liouville and other models. In that reference the first-order differential 
equations obeyed by the BPS solutions of the models are used to construct a mapping 
between the fields of the two models whose equivalence is to be demonstrated. In fact, one 
deforms one known model by using the mapping function and obtains another known model. 
It is shown that both models possess the same BPS energy when the deformation is performed 
in the Lagrangian density, in contrast to what happens when the deformation is carried out in 
the differential equations [13j. One can see that the BPS energy found in equation (IT^ of the 
model given by the Lagrangian density is the same of the model described by ([T]). Since 
we are not interested in proving the equivalence between those models we would rather prefer 
to work with an effective Lagrangian density whose structure is of the type kinetic— potential 
which leads to static Euler-Lagrangian equation of the type 0" = dUeff/d(j), whose classical 
solutions are also solutions of the equation ffTTl) . This procedure is more like the one of [13] 
and it is more convenient when one works with scalar fields in interaction with gravitation 
as we consider in the next section. 

2.1 Particular cases: 

Here we consider the same model with two interacting scalar fields applied in brane worlds 
scenario in the references [5] and [6] whose superpotential is 

w{^,x) = mY-^') + ^'^']■ ^^^^ 

In a recent paper, Dutra [lOj was able to show that for this case the orbit equations 
relating both fields can be obtained explicitly. They can be written as 

p{x) =<P'-a' = c^x^l^ - Y^X\ for A ^ 2/. (14) 

A — //X 
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and 

p(x) = 02 - = x^[Hx) + ci], for A = 2/i, (15) 

where Cq and Ci are constants of integration. In general, only the first of the above orbit 
equations is used to find well behaved classical solutions. One can check, for example, that 
the second of the above orbit equation fails to reproduce the minima of the model, namely 
= ±a and x = . From now one we consider a > 0. 

The first orbit equation can be used to construct models with only one scalar field which 
exhibit the main features of this model with two interacting scalar fields. An one example 
we consider the situation in which fi = X. In this case the orbit equation is given by 

+ CoX - (0' - a') = 0, (16) 

and the field x can be expressed in terms of as 



X = /(0) = -|±2V^o + 4(0^-«^)- (17) 

By substituting the above expression with the upper sign in the expression for PF0(0, x = 
/(0)) we obtain the following superpotential 

= 2/i(0' + - a' + + b^- (18) 

From now on we can work with a model described by the Lagrangian density 

£= V09^0-[/e//(0), (19) 

with the effective potential Ueff{<f)) = W^{(f))/2. In this model the first-order differential 
equations satisfied by 0(x) is 0' = and the solutions for the field are the same 

obtained in the model with two interacting scalar fields [6], [10] described by the Lagrangian 
density ([H) with W{(f),x) given in ( fT3l) . 

The choice of the lower sign in the expression ( |T71) would result in an effective potential 
Ueff{(j)) = with only one minimum and we are interested in effective potentials with 

at least two minima. The constant b = cq/2 must satisfy the inequality b < —a, such 
that we have non-singular solutions and the effective potential presents two global minima 
and one local minimum for a certain range of the parameter b. On the other hand, for 
b = —a, we have an effective potential with three global minima. We show in the Figure 1 
the behavior of the effective potential, in units of function of 0, for two different 

values of b and a = 1. One can see that for b < —a, the effective potential tends to exhibit 
three minima as b gets closer to —a. In fact, for b = —a the effective potential becomes 
Ueff{(f)) = 2/i^(0^ — a|0|)^that has also a minimum at = 0. The classical solutions of the 
first-order differential equation flTUl) . for b < —a, are given by 

^ ^ sinh(2/iax) 

cosh(2/xaa;) — b/ f ' 
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where / = ^W- — and the upper (lower) sign stands for the kink (anti-kink) solution. In 
the Figure 2 it is shown a profile of the topological classical solution for a sufficiently large 
value of and a double kink profile, usually called two-kinks, for h close to the critical 
value —a. This kind of configuration also arises as solution of the models introduced in the 
references [7] and [8]. 

Another choice of the parameter A, namely A = 4/i, leads to a similar non-polynomial 
effective potential. In this case, the orbit equation can be written as 



2 = 2. 

^ 4cr 



1 ± x/l + 16co(02 - a2)l , (21) 



and, by taking the upper sign in the above equation, we have the following effective super- 
potential 

= 4/i (^02 _ ^2 ^ rfcV0^ - d"^) , (22) 

where d = a/^/c^ — 1 and c = — IGcqci^. In order to have a well defined model, the 

constant cq must satisfy the inequality cq < l/16a^ and the solutions are given by 

^ ^ sinh(4/iax) ^^3^ 
cosh(4/iax) + c 

It worth mentioning that the model we have constructed in this case (A = 4/i) is very 
similar to the one proposed in the reference [7] if we set /i = 1/4 and make some identification 
between the parameter d with the parameter Bq from that paper. The behavior of the effective 
potential as a function of (j) is almost identical to the one presented in the previous case. 
Moreover, if cq = l/lGa^ one obtains the effective potential, Ueff{(j)) = 8/i^(0^ — a\(f)\)'^, with 
three minima. 

Effective polynomial potentials can also be obtained from the model described by the 
Lagrangian density ([T]) with W{(f), x) given by ffT^ . This is done by using the orbit equation 
(IT^ to express the field as a function of the field x- By conveniently rewriting the orbit 
equation as 

= gix) = ±^Cox'^^^-Y^x' + a^ (24) 
we have the effective potential given by 

UefAx) = \wl = 2/i\^ (co//'^ - + a') . (25) 



This potential has at least one minimum at x = 0, and we have to set X/fi = n, where n is 
a positive integer {n ^ 2), in order to have well defined potentials with, at least, two minima. 
This last condition leads to some constraints over the values of cq. We note that for A//i = 1, 
a positive definite potential with two minima, is obtained if cq = ±2a and, for A//i = 4, we 
find a positive definite potential with three minima at x = 0, ±2a if Cq = l/16a^. Those 
specific values for cq are the critical ones that lead to polynomials models with respect to the 
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field which possess three minima. We look for the critical values of cq, corresponding to 
different values of n, in order to have positive definite effective potentials, Ueff{x), with more 
than one minimum. First, we note that effective potentials of the type Ueff{(j)) = (0^ — a|0|)^ 
has classical solutions that connect the minimum = to the minima = +a or = —a 
and vice-versa. Moreover, we recall that those polynomial potentials, are constructed by 
substituting the orbit equation x = into the superpotential = A(0^ — a^) + fix^ and 
this one, by its turn, is substituted in Ueff{(p) = which is positive definite everywhere. 

If = is a minimum of this effective potential, it must correspond to x = ^^/na, due 
to the orbit equation. By substituting one of those values of x in the potential Ueff{x) 
and imposing that they are minima of this last, a priori, positive definite effective potential 
{Ueffi X = ±v^a) = 0), the critical value of cq 



Cq = a 



n~2 



n"/^ (26) 



is obtained. Except for the case n = 1 that presents two minima, one finds positive definite 
effective potentials with three minima only for n even. Such kinds of polynomial potentials 
were already discussed extensively in the literature |15]. These polynomial potentials have 
typical kink solutions that connect the minimum x = to one of the other two minima 
X = '^\fna and vice- versa. For n = 4, we have for instance 

cosh(/iax) ± sinh(/iax) 
A/cosh(2/iax) 

which have a kink profile very similar to that of the first-order differential equation solutions 
for the effective models of the kind ?7(0) = (0^ — ct|0|)^- These solutions that can be seen 
as half-torsion in a spin chain, also have similar profile to those exhibited by self-consistent 
solutions for inhomogeneous chiral condensates in the Nambu-Jona-Lasinio model in 1+1 
space-time dimensions [16]. Since these solutions are not continuously deformable into two- 
kinks solutions, they will not be applied to the brane-world scenario considered below where 
we are going to discuss the consequences of thicker branes in the warping of the space. 



3 Application to a brane world scenario 

We now consider the scalar field coupled to gravity in (4,1) space-time dimensions described 
by the action 

S = ^ d'xdy^\ [-\r + ^a„0a> - r(0)^ , (28) 

where g = Det{gab) and the metric is 

ds^ = gabdx^dx^ = e^^^'^^r^^^c/x^rfx^ - dr^, a,b = 0, 4, (29) 

where r = is the extra dimension, i]^^ is the Minkowski metric and e^^^^^ is the so-called 
warp factor, which is supposed to depend only on the extra dimension. The Greek indices 
run from to 3. 
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The static equations of motion following from the action (128|) and for the case that the 
scalar field depends only on the extra dimension are written as 



(ir^ 3 \ dr 

and 



dr'^ dr dr dcp 
d^A 2 fd(i)V 

[61) 



We consider that the potential V{(t>) can be written as [17] 



In this case the BPS solutions [IT] of the following first-order differential equations 

rf0 dWU) , dA 2^^,, , 
- ^ and — = T:^W{(f)) (34) 



dr d(j) dr 3 

are also solutions of the second-order differential equations f l30l) and fl3T]) and the equation 
( !32l) is identically satisfied. By taking VK</,(0) given by equation ( fTSj) together with the 
corresponding kink solution in ( l20l) . the superpotential is given by 



W{<p) = 2/i 



(35) 



and the warp factor is found by integrating the second of the equations ( IMl) with the classical 
solution (!20l) substituted in W{(j)). We show in Figure 3 two profiles of the warp factor 
(a = 1), corresponding to two different values of the parameter b; one close to and the other 
one far from the critical value b = —1. One can note that as far as b decreases (its modulus 
increases) the warp factor becomes more narrow. For values of b close to the critical value 
one can observe that the warp factor becomes wide and one can see a flat region, which 
signalizes a Minkowskian metric inside the domain wall. Thick is the brane, wider is the 
warp factor. We have also analyzed the case in which W{(f)) given by equation fl22l) and 
have verified that the behavior of the warp factor is not substantially different from that 
presented here. 

Now we consider the stability of the system by analyzing the equations of motion of linear 
small fluctuations around the classical solutions. This issue is also important to realize the 
localization of the gravity inside the domain wall [2]-[l7j. This is done by means of a 
perturbation of the metric, ds"^ = e'^^^'^\ri^y + eh^y)dx^dx^ — (ir^and a small perturbation 
around the classical solution — 0(r) +£0(r, x^), where e is a small number. By performing 
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those perturbations in the Lagrangian density and by expanding it up to 0{e), we obtain 
the following equations of motion 



dr dr d(j)^ 2 dr dr ^ 
for the fluctuation of the scalar field and 

+\ri,.e'^'^-^d^{t'K,) + = , (36) 

for the fiuctuations of the metric. 

In general it is quite difficult to take into account linear fiuctuations of all components 
of the metric together with the fiuctuations of the brane in order to have a broad view of 
the linear stability of the whole system. This is due to fact that the above set of coupled 
differential equations involving the fluctuations is very intricate to be solved. Nevertheless, 
it is possible to show that the transverse and traceless part of the fluctuations of the metric 
( h^y) decouple from fluctuations of the brane [T7]. By constructing h^^ = P^uaph"'^ from 
the projector operator P^^^p = \{'T(^aT^up + h^^uTTuo) - \T^iiuT^ap with vr^^ = r]^^ - d^d^/U, we 
have that the equation ( l36l) for the transverse and traceless part of the fluctuations of the 
metric simplifles to 

^ + - e-2^n V = 0. (37) 

dr^ dr dr 

By using separation of variables and expressing h^j, conveniently as 

V = e^'^''-'^e-(='/^)-^«e,.(r), (38) 

we flnd, by making use of the transformation of variables z = J e~^^^^dr, that ^fj,uiz) satisfles 
the following Schrodinger-like stability equation 

-%^ + V.ff{z)^,u = k%., (39) 

where 

, , 9 fdAV 3d^A 

is the effective potential. We show in the Figure 4 this effective potential against the variable 
r for two different values of b, one close and the other one far from the critical value. One 
can note that the shape of the effective potential for values of b far from the critical value is 
similar to the shape of others presented in the literature, for example in [2j. Therefore, the 
only bound-state solution is the one associated to the zero-mode which can be seen as the 
localization of the gravity inside the domain wall. Higher energy modes are non-localized 
states that can escape from the domain wall and propagate along the extra dimension. In 
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any case, no matter how big the value of |6| is, there is no room to have locahzed tachyon 
modes {k"^ < 0). The differential operator in equation (139|) can be factorized as the product 
of two operators which are adjoint of each other: 



3d A 



d 

dz 2 dz 



d 3dA\ 
'Tz ^ 2~dl) ^' 



fJ,U 



(41) 

|a|n)|^ > 0. Par- 

g(3/2)A(.) ^Yie 

corresponding transverse and traceless part of the fluctuation of the metric presents no de- 
pendence with the extra dimension. 

Another aspect that the non-polynomial models share with other models is the behavior 
of the matter energy density. It is given by 



then, for the n-th normalized eigenmode \n) we have fc^ = {n\a^a\n) 
ticularly, the non-normalized zero-energy eigenmode is given by 



(0), 



eir 



.4A{r) 



(42) 



where V{(j)) is the potential in ( 133|) evaluated at the classical solution. The behavior of the 
matter energy density is shown in Figure 5 for two different values of b. The features shared 
with most of the models is the peak of the energy density around the thick brane, which can 
be observed for values of b far from the critical value, and the presence of regions outside 
the domain wall where the matter energy density is negative. For values of b close to the 
critical one, two relative small peaks of the energy density show up around each wall of the 
double- wall structure, this is also a common feature of the models with thicker branes [EllH]- 
Although the energy matter is negative, the energy functional [T8] 



1 r+oo 

F= -l dr e^^^'-) 



(43) 



is positive definite. This energy functional generates the Euler-Lagrange equations (!30l) - (!32l) 
and is minimized by the solutions of the first-order differential equations (l34l) . since it can 
be rewritten as 



dr 



o4A(r) 



dr 



d 



±e^^(")W^(0)| 



(44) 



Note that the W^|0(oo) — W^|(/i(-oo) is bigger (less) than zero if we choose the upper (lower) 
sign. This energy plays the role of the topological Ebps energy in the present scenario of 
scalar fields interacting with gravitation in a (4,1) dimensional warped space-time. 

4 Final Remarks 



We have constructed effective models with only one scalar field which supports deformed 
kink solutions which has been called two-kink solutions. This has been done by means of 
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a general orbit equation relating two scalar fields and by eliminating one of them in terms 
of the other one. When such a procedure is done in the Lagrangian density, the effective 
Lagrangian density is not given in the usual kinetic— potential form (see equations ([S])-®)- 
In order to have an effective canonical Lagrangian density we adopt the elimination of one of 
the field in the equations of motion, instead. The model we have constructed here together 
with their two-kinks solutions are very similar to the one propose in [7], but, in principle, 
other (polynomial) models are obtained from that one we have started with. As a matter of 
fact, a myriad of models with two scalar interacting fields have been studied in reference [12] 
and from those we can, in principle, construct other effective models with only one scalar 
field. One of the models which is under analysis is the fourth model proposed in |19], whose 
potential is given by V{(f), x) = -^0 + fiX ~ '^/4(0^ + + 60^x^)- 

One of the effective models constructed here is applied to a brane-world scenario where 
its influence in the warp factor leads to a flat geometry inside the thicker domain walls. 
That model can also be used to analyze universal aspects of thick branes splitting in a 
warped geometry |3]. That phenomenon can be interpreted as a phase transition due to 
the variation of the parameters of the potential, which in our case are the parameters b 
for the superpotential f|T8l) . and c for the superpotential fl22l) . In each case there will be a 
transition in the form of the potential, from one with two vacua to another one with three 
vacua when the parameter reaches a critical value. One can think of b, for instance, as 
dependent on the temperature. Far from the critical temperature we have a single relatively 
thick brane and as the critical temperature is approached, the parameter b approaches —a 
from the left. For values of b close to but less than —a the thick brane splitting starts to 
take place (note the behavior of the solid curve in Figure 2) and the local minimum in the 
effective potential Ugff presents the tendency to become a global minimum. The splitting 
also influences the localization of the zero modes fluctuation of the metric as can be seen 
by the behavior of the effective potential in Figure 5. In that flgure one can see that for 
values of b far from the critical value the effective potential has a volcano shape which 
leads to a narrow (more localized) zero-mode fluctuation than in the situation in which b 
is close to the critical value. The splitting of the brane is also manifested in the splitting 
of the matter energy density shown in Figure 6. When b = —a, the effective potential is 
Ueff{(p) = 2/i^(0^ — a|0|)^ and the brane is split in two. We can suppose that the branes are 
at a distance 2L from each other and localized around the core of each one of the solutions 
0_(r) = -a/2{l - tanh[2/ia(r + L)]} and 0+(r) = a/2{l + tanh[2/xa(r - L)]. We have 
not analyzed what happens to the warp factor and the fluctuation of the metric under the 
influence of both scalar flelds simultaneously but this might be done by following the same 
numerical approach adopted in the reference [2U], where it was obtained the spectrum of 
fermions in the background of a kink and a anti-kink which are far apart from each other. 
We think that the effective potential for the fluctuations of the metric will have the shape 
of two volcanos whose craters are distant 2L from each other. That effective potential can 
support a zero eigenmode and the eigenfunctions may have peaks around each brane or in 
the region between the branes. 

The effective models presented here can also be considered in space-time with D > 4 
dimensions and with the solutions depending only on the radial coordinate. Since nonlinear 
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models with only scalar fields are not stable in space-time dimensions bigger than two, as 
has been demonstrated by Derricks theorem [21], one should resort to a convenient bypass 
by introducing an explicit dependence of the interacting potential on the coordinates, as has 
been provided by [5]. The resulting radial solutions are very similar to those shown in the 
reference [7]. Those solutions and their consequences in warped space-time with two and 
three extra dimensions is under study and the results will be reported elsewhere. 
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Figure 1: The effective potential in tlie case A = /x, for a = 1 and b = —1.001 (solid line) 
and b — —1.3 (dashed line). 
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Figure 2: Typical kink-profilcs, for a = 1 and b — —1.001 (solid line), corresponding to 
two-kinks solution and b — —1.3 (dashed line), corresponding to a single kink solution. 
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Figure 3: Warp factor for two different values of b. One closed to the critical value (thick 
line) with a meseta shape and the other one far from the critical (thin line). 
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Figure 4: The effective potential in the effective Schrodinger equation for the fluctuation of 
the metric for two different values b, b — —1.0001 (thick line) and b — —1.2 (thin line). 
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Figure 5: The matter energy density in the case X — fx, for a — 1 and b — —1.0001 (thin 
hne) and b — —1.3 (thick hne). 
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